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ABSTRACT 


This  thesis  deals  with  the  characterization  of  the  dependence  on  the  flow 
geometry  of  the  streamwise  fluctuations  of  the  stagnation  point  of  vortex  breakdown  in 
axisymmetric  tubes  and  over  delta  wing  aircraft.  The  statistical  analysis  presented  herein 
shows  that  in  an  axisymmetric  tube  the  ’darting’  about  the  mean  stagnation  point  are 
distributed  normally  for  the  two  Reynolds  numbers:  Reo  =  230,000  and  300,000 
(independently  of  the  Reynolds  number  in  the  range  noted).  The  darting  over  a  delta 
wing  is  not  only  non-Gaussian  but  also  exhibits  rather  large  localized  fluctuations 
(Strouhal  numbers  ranging  from  0.04  to  0.1),  presumably  due  to  the  strong  influence  of 
the  surrounding  flow  and  the  geometrical  conditions:  increase  of  circulation  along  the 
trailing  edge,  the  abrupt  separation  of  flow  at  the  base  of  the  delta  wing,  and  other 
protuberances  that  emerge  from  the  upper  and  lower  surfaces  of  the  wing  (support 
elements  in  laboratory  and  stabilizers  on  delta  wing  aircraft).  It  is  concluded  that  the 
behavior  of  vortex  breakdown  is  strongly  dependent  on  the  surrounding  geometry  and 
that  only  experiments  in  axisymmetric  tubes  can  provide  the  purest  form  of  vortex 
breakdown  for  numerical  simulations  and  analytical  studies  towards  the  understanding  of 
the  internal  turbulence  and  its  spectrum  within  the  breakdown  bubble  for  theoretical  and 
industrial  purposes. 
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I. 


VORTEX  BREAKDOWN 


A.  INTRODUCTION 

1.  Background 

Vortex  breakdown  is  a  phenomenon  that  has  been  studied  extensively  both 
theoretically  and  experimentally,  yet  it  remains  one  of  the  many  such  naturally  occurring 
phenomena  in  fluid  flows  that  has  not  yet  been  fully  understood.  As  is  the  case  of  many 
other  flows,  researchers  have  not  been  able  to  acquire  sufficient  understanding  of  vortex 
breakdown  over  delta  wings,  in  tubes,  and  in  combustion  chambers.  Observations  have 
been  made  that  describe  the  fundamental  behaviors  of  vortex  breakdown  using  many 
different  and  varying  experimental  procedures  and  equipment.  But  prior  to  exploring 
these  rather  complex  methods  for  categorizing  vortex  breakdown,  it  is  important  to 
acquire  a  rather  basic  understanding  of  the  physics  and  forms  of  the  flow  through  an 
investigation  of  the  breakdown  in  tubes. 

The  most  common  tool  for  identifying  structures  in  fluid  flow  is  called  flow 
visualization.  Flow  visualization  encompasses  introducing  a  visually  contrasting  agent, 
dye,  into  the  flow  in  order  to  trace  streamlines  within  the  flow.  In  vortex  breakdown 
(hence  forth  VB)  the  dye  is  injected  through  a  port  aligned  with  the  axis  of  the  vortex 
core.  As  the  fluid  progresses  naturally  toward  the  breakdown  point,  the  dye  defines  the 
flow  structures  in  the  fluid.  This  kind  of  flow  visualization  enables  one  to  describe 
various  aspects  of  the  flow  that  evade  state  of  the  art  instrumentation  such  as  laser 
Doppler  velocimetry  (LDV)  and  digital  particle  image  velocimetry  (DPIV)  which  are 
used  to  measure  quantitative  elements  of  the  flow.  Other  such  instrumentation  is  used  for 
detennining  the  real-time  behavior  of  the  flow  as  it  applies  to  aspects  of  mixing  in 
combustion  chambers  and  lift  and  drag  in  aerodynamics. 

The  two  fundamental  parameters,  in  addition  to  geometrical  variables,  that  govern 
the  behavior  of  the  flow  are  Reynolds  and  swirl  numbers,  defined  by  equations  (1)  and 
(2)  respectively. 


(1) 


Re0 


UqDq 

v 


1 


r 


(2) 


r  = 


u0d0 


Where  Uo  is  the  mean  axial  velocity,  Do  is  the  tube  entrance  diameter,  v  is  the  kinematic 
viscosity,  and  T  is  the  circulation. 

Because  of  the  dependence  of  these  parameters  on  geometry,  they  are  oftentimes 
presented  in  different  forms.  As  such,  they  serve  as  scaling  parameters  and  can  be  used 
to  compare  data  with  those  obtained  by  others.  The  current  investigation  employs  a  tube- 
and-vane-type  apparatus  that  allows  for  the  precise  setting  of  the  desired  Reynolds  and 
swirl  numbers  through  the  regulation  of  the  flow  rate  and  the  adjustable  vane  angles.  The 
combination  of  the  flow  visualization  with  the  variability  of  these  experimental 
conditions  provides  insight  into  what  happens  as  VB  occurs. 


2.  Forms  of  Vortex  Breakdown 

Understanding  the  forms  of  VB  is  important  in  determining  specific 
characteristics  and  behaviors  of  the  flow.  It  is  also  important  to  note  the  relationship 
between  the  fonn  of  the  breakdown  and  the  conditions  that  are  present  for  that 
breakdown  to  occur,  i.e.  Reynolds  and  swirl  numbers.  The  four  major  types  of  VB  are 
the  double  helix,  spiral,  nearly  axisymmetric  recirculation,  and  conical  breakdowns. 
Several  combinations  of  these  salient  forms  of  breakdown  also  transpire  to  create  a 
medley  of  structures  interacting  within  the  flow. 

a.  Helical  Vortex  Breakdown 

Understanding  the  major  forms  of  VB  demands  an  exploration  of  the 
literature  that  helped  establish  them.  The  double  helix  was  first  observed  by  Sarpkaya 
(1971a).  He  used  a  tube-and-vane-type  apparatus  with  adjustable  flow  rate  and  vane 
angles  to  control  the  Reynolds  number  and  swirl  number.  Dye  injected  along  the 
centerline  into  the  flow  was  observed  to  decelerate  and  expand  into  two  thin  triangular 
sheets  that  wrap  around  each  other  forming  a  double  helix.  This  form  gradually  expands 
to  the  limits  of  the  tube  and  breaks  down  into  mild  turbulence.  Sarpkaya  (1971b)  notes 
that  the  double  helix  is  highly  unstable  and  reacts  to  internal  and  external  disturbances. 
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For  this  reason  Sarpkaya  (1971b)  states  that  the  double  helix  may  not  be  a  form  of  vortex 
breakdown,  but  rather  a  “consequence  of  the  instability  of  the  prevailing  flow  to  spiral 
disturbances.” 

Figure  1  illustrates  this  fonn  of  VB.  The  double  helix  is  uniquely  different 
from  the  other  forms  of  VB  because  the  flow  remains  laminar  throughout  the  breakdown. 
Figure  1  below  shows  the  dye  filament  at  the  vortex  core  splitting  into  two  triangular 
sheets  that  become  intertwined  upstream  of  the  shearing  point.  The  sheets  expand  and 
become  thinner,  making  them  susceptible  to  turbulent  motions.  It  is  only  near  the  walls 
of  the  diverging  tube  that  sheets  become  interrupted  by  turbulence. 


Figure  1.  Double  Helix  vortex  breakdown  in  a  diverging  tube.  (After  Sarpkaya 

1971a) 

The  direction  of  the  swirling  flow  is  a  consequence  of  the  inception  of  the  instability  and 
of  little  importance  to  the  breakdown.  Both  sheets  curl  in  towards  one  another  regardless 
of  the  direction  of  swirl.  In  the  top  figure  of  Figure  1  the  sheets  curl  out  of  the  page, 
while  in  the  bottom  figure  the  double  helix  is  curling  in  the  opposite  direction.  Because 
of  the  relatively  low  Reynolds  numbers  and  high  circulation  numbers,  reportedly  Re  < 
2000  and  Cl  >  2.3  from  Sarpkaya  (1971a),  the  double  helix  has  only  been  observed  in 
laboratory  settings. 
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b.  Spiral  Vortex  Breakdown 

Spiral  vortex  breakdown  is  the  next  prominent  type  of  vortex  breakdown, 
occurring  at  Reynolds  and  swirl  numbers  slightly  greater  than  those  for  the  double  helix 
form.  Sarpkaya  (1971a)  reports  that  at  constant  Reynolds  number  of  about  5,000  the 
vane  angles  can  be  increased  to  about  20  degrees  to  see  the  first  indication  of  VB.  By 
increasing  the  vane  angle  farther,  the  vortex  core  dye  filament  breakdown  moves  rapidly 
upstream  and  develops  an  abrupt  kink  that  is  indicative  of  rapid  deceleration  of  the  flow 
at  the  vortex  core.  After  one  or  two  windings  the  filament  breaks  up  into  large-scale 
turbulence.  Figure  2  shows  the  spiral  breakdown  observed  by  Sarpkaya  (1971a). 


Figure  2.  Spiral  vortex  breakdown  in  a  diverging  tube.  (After  Sarpkaya  1971a) 


Figure  2  shows,  from  top  to  bottom,  the  effect  of  increasing  the  swirl.  In  each  case,  the 
nascent  spiral  is  followed  by  large-scale  turbulence.  More  detailed  observations  show 
that  turbulent  mixing  increases  along  the  length  of  the  tube.  It  is  also  apparent  from 
Figure  2  that  the  kink  in  the  dye  filament  is  not  stationary.  This  unsteadiness  is  one  of 
the  behaviors  of  VB  that  is  yet  to  be  quantified.  The  changes  in  location  of  what  may  be 
called  the  stagnation  point  make  the  accurate  quantification  of  time-averaged  statistics 
rather  difficult. 
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c.  Nearly  Axisymmetric  Recirculation  Vortex  Breakdown 
The  characteristic  shape  of  the  axisymmetric  recirculation  form  of  VB  is  a 
'bubble'.  This  description  may  be  misleading  to  the  casual  reader,  since  shape  and 
function  of  the  shape  are  two  separate  entities.  Therefore,  it  is  best  to  introduce  this  form 
of  VB  with  the  help  of  a  figure.  Figure  3  shows  what  happens  when  the  swirl  is 
sufficiently  increased.  When  the  dye  filament  in  the  vortex  core  reaches  the  stagnation 
point  it  spreads  around  the  outer  shell  of  the  bubble.  Near  the  downstream  limit  of  the 
shape,  the  dye  is  drawn  into  the  recirculation  portion  of  the  bubble  and  becomes 
entrained  there,  as  shown  in  the  bottom  figure  of  Figure  3. 


Figure  3.  Axisymmetric  recirculation  “bubble”  vortex  breakdown.  (After  Sarpkaya 

1971a) 

It  comes  as  no  surprise  that  this  kind  of  recirculation  could  be  utilized  in  combustion 
processes  to  improve  efficiency  and  to  reduce  combustion  products.  Lilley  (1977)  and 
Swithenbank  and  Chigier  (1969)  provide  comprehensive  reviews  of  the  experimental 
work  done  in  this  field. 

Unlike  spiral  VB,  axisymmetric  recirculation  VB  does  not  initiate  large- 
scale  turbulence  in  the  adjacent  regions.  Figure  4  shows  a  sheet  of  dye  that  surrounds  the 
bubble,  expanding  to  accommodate  the  displaced  volume  of  the  bubble  and  converging 
again  after  the  bubble. 
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Figure  4.  Axisymmetric  vortex  breakdown  with  a  laminar  sheet  surrounding  the 

vortex  breakdown.  (After  Sarpkaya  1971a) 

d.  Conical  Vortex  Breakdown 

In  1995,  a  fourth  form  of  VB  was  discovered  and  investigated  by 
Sarpkaya  (1995).  Since  then  the  conical  VB  has  been  the  focus  of  the  most  recent 
research  in  the  field.  Occurring  at  high  Reynolds  numbers  in  noncavitating,  turbulent, 
swirling  flows,  this  type  of  breakdown  exhibits  a  clear  stagnation  point  followed  by  an 
axisymmetric  expanding  cone.  Sarpkaya’ s  (1995)  observations  show  a  small  area  of 
recirculation  immediately  after  the  stagnation  point  that  bursts  into  a  conical  spiraling 
wake.  By  increasing  the  Reynolds  number  based  on  the  tube  diameter  to  200,000  and 
adjusting  the  circulation  number,  he  was  able  to  eliminate  the  axisymmetric  recirculation 
before  the  cone.  Figure  5  shows  the  progression  of  the  VB  with  increasing  Reynolds 
number.  The  Reynolds  numbers  are  50  000,  100  000,  150  000,  and  200  000  for  the  top  to 
bottom  pictures  respectively.  The  corresponding  circulation  numbers  are  0.77,  0.61, 
0.55,  and  0.50  from  the  top  to  the  bottom  picture,  respectively. 
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Figure  5. 


Conical  vortex  breakdown.  (After  Sarpkaya  1995.) 


Novak  and  Sarpkaya  (2000)  discuss  two  phenomena  that  occur  with  the  conical  VB, 
vortex  core  meandering  and  stagnation  point  darting.  The  first,  vortex  core  meandering, 
occurs  when  the  core  visibly  moves  radially  about  the  mean  centerline.  Novak  and 
Sarpkaya  (2000)  state  that,  when  collecting  data  using  an  LDV,  core  meandering  creates 
a  region  where  turbulence  measurements  become  unreliable.  The  darting  of  the 
stagnation  point  of  the  conical  VB  is  the  second  phenomena  and  will  encompass  the  bulk 
of  the  discussion  in  this  thesis.  It  is  a  random  streamwise  fluctuation  of  the  location  of 
the  VB,  with  no  noticeable  transition  time  or  periodic  motion.  There  is  however  a  mean 
location  for  the  stagnation  point  and  a  normal  distribution  of  the  axial  displacements  from 
the  mean,  as  shown  later. 

3.  Scope  of  the  Present  Work 

So  far  the  discussion  has  focused  primarily  on  experimental  flow  visualization 
results  from  tube-and-vane  type  apparatus  to  describe  VB.  Vast  amounts  of  data  exist 
that  examine  time-averaged  flows  and  provide  conclusions  that  relate  experimental 
results  to  applicable  real  life  scenarios.  By  utilizing  such  results,  engineers  are  able  to 
design  control  systems  in  platforms  such  as  the  delta  wing  to  provide  for  safer  operation 
and  enhanced  maneuverability.  The  present  study  is  aimed  at  furthering  the  bank  of 
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knowledge  that  engineers  have  to  work  with.  By  understanding  the  behavior  of  the 
darting  of  the  stagnation  point  of  VB,  engineers  should  be  able  to  take  one  step  closer  to 
developing  a  safe  operating  envelop  for  the  delta  wing  aircraft. 


B.  PHENOMENA 

Experiments  dealing  with  the  fluctuations  of  VB  over  delta  wings  have  shown 
that  a  large  number  of  parameters  must  be  taken  into  consideration  to  understand  the 
sources  as  well  as  the  physics  of  the  phenomenon.  These  include,  but  are  not  limited  to, 
the  sweep  angle,  leading-edge  geometry,  wing  thickness,  angle  of  attack,  and  the  free- 
stream  conditions.  On  the  other  hand,  the  fluctuations  of  VB  in  tubes  are  controlled  by 
fewer  parameters  such  as  the  Reynolds  number  and  the  circulation  of  the  upstream 
vortex.  There  are,  to  be  sure,  secondary  influences  such  as  the  boundary  layer  growth  on 
the  tube  wall  and  the  effect  of  air  bubbles  for  the  experiments  conducted  using  water  as 
the  working  fluid.  Nevertheless,  the  tube  experiments  provide  a  clearer  understanding  of 
the  characteristics  of  VB  in  general  and  of  the  streamwise  fluctuations  in  particular. 
Furthermore,  the  effect  of  the  boundary  layer  and  the  dissolved  air  can  easily  be 
minimized  by  the  slight  expansion  of  the  tube  wall,  by  de-aerating  the  water  before 
introducing  it  into  the  test  chamber,  and  by  conducting  the  experiments  under  sufficiently 
high  pressure,  as  in  the  case  of  the  experiments  reported  herein.  These  will  be  discussed 
in  greater  detail  later  following  a  brief  summary  of  the  delta  wing  experiments. 

For  delta  wings,  the  sweep  angle  and  angle  of  attack  are  two  of  the  parameters 
that  determine  where  the  VB  occurs.  Sweep  angle  is  a  measurement  of  the  angle  the 
leading  edge  makes  with  a  line  perpendicular  to  the  axis  of  the  delta  wing  that  passes 
through  the  apex.  Aspect  ratio  (the  ratio  of  the  wing  span  squared  to  the  area  of  the 
wing)  is  another  way  to  represent  the  sweep  angle.  Figure  6(a)  is  provided  for 
clarification  of  the  sweep  angle.  Figure  6(b)  shows  the  angle  of  attack  as  the  angle 
measured  from  the  plane  in  the  direction  of  the  freestream  fluid  motion,  usually  the 
horizontal  plane,  to  the  bottom  surface  of  the  delta  wing. 
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Figure  6.  (a)  Top  view  rendering  of  a  delta  wing.  The  sweep  angle  is  labeled. 

(After  Hummel  1967)  (b)  Side  view  rendering  of  a  delta  wing  showing  angle  of 

attack.  (After  LeMay  1990) 

Mitchell  (2000)  observed  that  “increasing  angle  of  attack  adds  energy  to  the 
leading-edge  vortices  as  a  result  of  increased  adverse  pressure  gradients  producing  higher 
axial  and  swirl  velocities  than  at  lower  incidence  angles.”  The  breakdown  responds  to 
these  increased  axial  and  swirl  velocities  by  moving  upstream  towards  the  apex  as  the 
angle  of  attack  is  increased.  Sarpkaya  (1971a)  observed  similar  events  in  a  tube  when 
increasing  the  vane  angle  to  increase  swirl.  In  his  tube-vane -type  apparatus,  the  increase 
in  swirl  resulted  in  an  axial  upstream  displacement  of  the  vortex  breakdown. 

Wentz  (1971)  found  that  increasing  sweep  angle  would  increase  the  angle  of 
attack  necessary  for  breakdown  to  occur.  At  a  threshold  of  75  degrees,  the  sweep  angle 
was  determined  to  have  no  effect  on  the  angle  of  attack  for  breakdown  to  occur.  Wentz 
(1971)  also  stated  that  no  breakdown  was  observed  for  a  sweep  angle  of  45  degrees  at 
any  angle  of  attack  and  at  sweep  angles  of  50  and  55  degrees  the  breakdown  never 
reached  the  trailing  edge.  This  puts  into  focus  that  delta  wings  with  sweep  angles  of 
greater  than  55  degrees  will  experience  VB  over  their  wings,  resulting  in  changes  to  lift, 
drag,  and  maneuverability.  Evidently,  the  frequency  and  amplitude  of  the  fluctuations  of 
the  location  of  the  breakdown  can  be  critical  for  stability  of  delta  wing  aircraft. 

One  critical  angle  of  attack  will  be  the  angle  of  attack  where  the  VB  occurs  over 
the  trailing  edge  of  the  aircraft.  Hummel  (1967)  and  Mitchell  (2000)  discuss  the 
conditions  for  the  VB  to  occur  over  the  delta  wing  in  detail.  Results  from  both  studies 
show  that  as  the  sweep  angle  is  increased,  the  angle  of  attack  for  VB  to  occur  over  the 
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wing  also  increases.  Mitchell’s  (2000)  results  also  show  that  as  the  freestream  velocity 
increases,  the  angle  of  attack  for  VB  to  occur  over  the  wing  increases.  Thus,  it  may  be 
concluded  that  an  increase  in  freestream  velocity  and/or  sweep  angle  will  increase  the 
necessary  angle  of  attack  for  VB  to  occur  over  delta  wings. 

Hummel’s  (1967)  investigation  also  includes  results  for  the  lift,  drag,  and  pitching 
moment  coefficients.  Once  the  VB  crosses  the  trailing  edge  and  is  over  the  wing,  there  is 
a  marked  decrease  to  these  coefficients.  As  the  breakdown  moves  toward  the  apex,  the 
condition  becomes  more  serious  and  the  delta  wing  begins  to  experience  stalling. 
Hummel  (1967)  concluded  that  “at  larger  aspect  ratios,  the  effect  of  vortex  breakdown 
starts  at  lower  values  of  the  angle  of  incidence.”  This  conclusion  is  complimentary  to  the 
one  made  previously  and  simply  states  that  as  the  sweep  angle  is  decreased  (aspect  ratio 
increases)  the  VB  will  cross  the  trailing  edge  at  a  smaller  angle  of  attack. 

Menke  (1996)  offered  further  validity  to  the  discussion  of  Reynolds  number 
dependence  by  noting  that  in  his  frequency  analysis  the  amplitude  of  the  dominant 
frequency  of  the  darting  of  the  stagnation  point  was  smaller  at  lower  Reynolds  numbers. 
This  study  also  investigated  the  amplitude  of  trailing  edge  pressure  fluctuations  for  three 
mean  breakdown  locations,  at  the  trailing  edge,  near  the  apex,  and  at  the  mid-point. 
Understanding  the  pressure  fluctuations  leads  to  an  understanding  of  the  stability  of  the 
aircraft.  For  the  breakdown  over  the  trailing  edge,  the  least  amount  of  fluctuation  was 
observed.  When  the  breakdown  was  over  the  mid-point,  pressure  fluctuations  at  the 
trailing  edge  increased  to  two  or  three  times  the  magnitude  of  those  observed  for 
breakdown  over  the  trailing  edge.  The  highest  fluctuation  of  pressure  at  the  trailing  edge 
was  observed  for  breakdown  close  to  the  apex.  These  fluctuations  were  two  to  three 
times  larger  than  those  for  the  mid-point.  As  these  fluctuations  increase  instability  of  and 
loss  of  control  of  the  aircraft  are  more  likely  to  occur.  Menke  (1996)  also  investigated 
the  distribution  of  the  fluctuations  of  the  breakdown  location  based  on  two  angles  of 
attack.  The  results  were  presented  as  a  histogram  and  a  probability  density  function. 
Menke  (1996)  concluded  that  his  results  were  not  normally  distributed  when  compared  to 
a  Gaussian  curve.  The  consequence  of  this  finding  is  that  the  fluctuation  in  VB  is  neither 
periodic  nor  normally  distributed.  This  point  will  be  taken  up  again  in  connection  with 
the  discussion  of  the  tube  experiments. 
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Menke  (1996)  also  investigated  the  interaction  of  the  vortices  occurring  on 
opposite  leading  edges.  He  stated  that  “there  is  an  antisymmetric  motion  of  breakdown 
location  for  left  and  right  vortices.”  Both  breakdowns  oscillate  at  about  the  same 
frequency  separated  by  a  phase  angle  of  approximately  180  degrees.  Mitchell  (2000) 
noted  that  there  are  strong  interactions  between  the  leading  edge  vortices,  such  that  “as 
the  breakdown  location  on  one  side  of  the  delta  wing  moves  aft,  the  other  tends  to  move 
forward  and  vice  versa.”  Both  Mitchell  (2000)  and  Menke  (1996)  are  in  agreement  that 
in  some  instances  the  oscillations  are  180  degrees  out  of  phase. 

The  brief  summary  here  illustrates  the  complexity  of  experimentation  with  models 
of  delta  wing  aircraft.  Quite  simply  it  is  the  shape  of  the  delta  wing  and  the  variability  of 
the  phenomena  that  makes  the  measurements  difficult.  What  is  even  more  difficult  is 
developing  a  clear  understanding  of  the  physics  of  the  fluid  motion.  That  is  why  the 
remainder  of  this  discussion  will  be  focused  on  breakdown  in  tubes. 

C.  CONICAL  VORTEX  BREAKDOWN 

1.  Characteristics  of  Conical  Vortex  Breakdown 

In  the  foregoing  we  have  discussed  the  fundamental  forms  of  VB.  We  have  stated 
that  where  it  occurs,  why  it  occurs,  and  what  happens  when  it  occurs  are  questions  that 
have  only  been  answered  partially.  Of  the  forms  discussed  only  the  conical  VB  (CVB) 
will  be  examined  here  to  uncover  answers  to  these  questions.  CVB  owes  its  occurrence 
to  high  Reynolds  number  but  not  necessarily  to  high  circulations.  The  fact  that  it  does 
not  occur  at  a  fixed  point  in  a  tube  nor  on  a  delta  wing  is  the  reason  for  the  work 
described  here  in.  Our  purpose  was  to  understand  this  aspect  of  CVB  to  the  extent 
possible  through  the  analysis  of  high  speed  motion  pictures  depicting  the  motion  of  CVB 
and  carrying  out  a  thorough  statistical  analysis  to  uncover  the  physics  of  the  phenomena. 
These  results  will  subsequently  be  compared  with  those  obtained  on  delta  wings  at 
various  Reynolds  numbers. 

As  shown  in  Figure  5,  CVB  starts  from  a  stagnation  point  and  expands  conically 
toward  the  limits  of  the  tube.  In  the  case  of  flow  over  a  delta  wing,  the  limits  are 
determined  by  geometry,  Reynolds  number,  and  angle  of  attack  to  name  a  few.  Because 
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of  the  nature  of  this  flow  it  is  no  longer  possible  to  identify  spatially  coherent  structures 
within  the  breakdown  region.  What  is  seen  are  only  various  scales  of  micro  to  macro 
turbulence  emerging  from  the  stagnation  point.  The  stagnation  point  itself  is  at  an  instant 
a  location  from  which  CVB  originates,  however  its  location  is  not  stagnant  in  the  strict 
sense  of  the  word.  Like  the  flow  inside  the  cone,  the  stagnation  point  location  is 
constantly  changing.  It  will  be  shown  later,  however,  that  there  is  order  to  the  chaos  in 
the  stagnation  point  location,  giving  physical  meaning  to  the  CVB  phenomena. 

2.  Streamwise  Oscillations  of  Conical  Vortex  Breakdown 

This  study  examines  CVB  at  two  Reynolds  numbers,  ReD=230  000  and  300  000, 
in  a  diverging  tube.  The  circulation  was  adjusted  such  that  the  mean  location  of  the 
breakdown  was  within  the  test  section,  but  not  to  the  point  for  cavitation  to  occur  in  the 
flow.  The  apparatus  was  also  pressurized  to  maintain  sufficient  static  pressure  to  prevent 
cavitation.  Aside  from  cavitation,  there  are  other  external  factors  that  require 
considerable  attention  and  are  worth  noting.  Because  the  flow  was  driven  by  a  high 
volume  pump,  the  variation  in  electrical  input  and  water  pressure  could  cause  noticeable 
deviation  in  results  if  that  variation  was  sufficiently  large.  Because  the  temporal  average 
of  velocities  and  turbulence  intensities  were  taken  over  relatively  long  periods,  small 
variations  in  electrical  input  and  water  pressure  would  not  present  a  notable  discrepancy 
in  the  results.  It  can  also  be  reasoned  that  in  the  case  of  larger  external  variations,  the 
response  of  the  VB  would  be  significant  enough  to  be  eliminated  on  the  basis  of  it  being 
an  outlier  in  the  data.  In  the  moderate  case  that  an  external  disturbance  was  promulgated 
to  the  results,  an  analysis  of  the  frequency  spectrum  would  help  separate  the  noise  from 
the  actual  behavior  of  the  VB.  Figure  7  shows  the  power  spectrum  density  for  the  darting 
of  the  stagnation  point  of  CVB. 
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(d) 

Figure  7.  Power  Spectrum  Density  for  (a)  450  data  points  at  ReD=230,000,  (b)  507 
data  points  at  ReD=300,000,  (c)  Non-dimensional  Power  Spectrum  Density  for 
Reo=230,000,  and  (d)  Non-dimensional  Power  Spectrum  Density  for 

Reo=3 00,000. 
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When  compared  to  results  in  the  literature  from  Menke  (1996)  and  Mitchell  (2000), 
whose  experiments  were  conducted  over  delta  wing  models,  the  peak  frequency  for  the 
oscillations  occur  at  the  same  non-dimensional  frequency  or  Strouhal  number  range,  from 
0.04  to  0.1,  as  seen  in  Figure  7(c)  and  (d).  The  absence  of  any  large  peaks  at  higher 
frequency  in  Figure  7(a)  and  (b),  such  as  60  FIz  for  electricity,  indicates  that  the  CVB 
fluctuation  is  free  from  such  external  disturbances. 


15 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


16 


II.  EXPERIMENTAL  EQUIPMENT  AND  PROCEDURES 


A.  FLOW  APPARATUS 

The  flow  apparatus  consisted  of  a  removable  Plexiglas  tube,  the  vane  section,  the 
foam  baffle,  steel  outer  chamber  with  four  windows,  inlet  and  outlet  ports,  and  a  concave 
centerbody  as  shown  in  Figure  8  (for  additional  details  see  Novak  1998). 
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Figure  8.  Tube-and-vane  type  flow  apparatus  for  inducing  vortex  breakdown  (From 

Novak  1998) 

Other  requisite  elements  were  the  centrifugal  pump,  air  escape  lines,  flow  meters, 
and  a  12  m  reservoir.  Each  of  these  parts  was  vital  in  maintaining  steady  conditions  for 
the  VB  to  occur.  Because  the  working  fluid  was  water,  eliminating  air  bubbles  and 
cavitation  became  a  design  consideration  when  working  at  high  Reynolds  numbers. 
Another  design  consideration  was  the  development  of  a  boundary  layer  on  the  inner 
surface  of  the  tube. 

Because  of  the  possibility  of  cavitation  and  for  air  bubbles  to  form,  the  steel  outer 
chamber  was  constructed  to  withstand  high  static  pressures.  Sufficient  static  pressure 
was  maintained  in  the  apparatus  to  prevent  cavitation  by  controlling  the  flow  rates  at  the 
inlet  and  outlet.  Additionally,  air  was  periodically  purged  from  the  system  to  reduce  the 
occurrence  of  bubbles  in  the  test  section. 
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Since  a  developing  boundary  layer  on  the  inner  surface  of  the  tube  was 
unavoidable,  the  tube  was  designed  to  account  for  the  displacement  thickness  of  the 
boundary  layer  by  diverging  at  a  1.4  degree  angle.  This  prevented  the  core  flow  from 
being  adversely  effected  by  the  developing  boundary  layer.  Figure  9  is  a  rendering  of  the 
tube  used  in  the  present  study. 


Figure  9.  Profile  of  test  tube  used  (all  dimensions  are  in  millimeters).  (After  Novak 

and  Sarpkaya  2000) 

More  specifically,  the  tube  depicted  in  Figure  9,  starts  at  the  inlet  radius  of  23 
mm,  diverges  to  a  radius  of  25.4  mm,  and  then  converges  to  a  throat  with  a  radius  of  23.3 
mm  located  at  140  mm  from  the  entrance.  From  this  point,  the  tube  rapidly  diverges  to  a 
28.5  mm  radius  at  223  mm  from  the  inlet  and  then  gently  diverges  at  1.4  degrees  to  the 
final  diameter  of  28.6  mm.  This  tube  was  machined  from  a  solid  three  inch  Plexiglas  rod. 
The  factory  finish  was  maintained  on  the  outer  surface,  but  the  inner  wall  was  machined 
to  the  desired  profile,  ground  to  eliminate  sharp  comers,  and  then  polished  to  a  mirror 
finish. 

Another  important  part  of  the  apparatus  is  the  concave  centerbody  which 
facilitated  the  transition  of  swirling  flow  from  the  vanes  into  the  test  section.  The 
centerbody  was  also  the  location  of  the  two  inlet  ports  for  the  flow  visualization  agents, 
food  coloring  and  fluorescein.  One  port  was  located  at  the  centerline  so  that  dye  could  be 
injected  into  the  vortex  core  and  the  other  port  was  located  slightly  off  the  flow  axis.  It 
should  be  apparent  that  the  port  along  the  centerline  was  used  to  obtain  flow  visualization 
of  the  VB,  which  was  of  paramount  importance,  but  the  other  port  provided  for 
observations  of  the  behavior  of  the  flow  surrounding  the  VB. 
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The  system  was  cleaned  periodically.  Before  a  new  run  commenced,  the  test  tube 
was  removed  to  facilitate  cleaning  of  the  inside  wall  and  the  vanes  and  bellmouth  were 
also  wiped  clean.  The  water  was  also  filtered  through  a  10  micron  filter  before  and  after 
the  experiments  to  remove  any  particulate  matter. 

The  apparatus  described  herein  was  an  augmentation  from  that  used  by  Sarpkaya 
(1971a,  1971b,  1974,  1995)  and  identical  to  that  used  by  Novak  and  Sarpkaya  (2000). 


B.  FLOW  VISUALIZATION  AND  OBSERVATIONS 

The  SVHS  video  that  was  analyzed  for  this  thesis  was  recorded  using  a  Sony 
Shoulder  mounted  camera  at  the  standard  speed  of  30  frames  per  second  for  CVB  with 
Reynolds  numbers  of  230,000  and  300,000.  Once  the  correct  flow  rates  were  achieved, 
the  vane  angles  were  adjusted  such  that  the  CVB  occurred  in  the  field  of  view  of  the 
camera.  The  amount  of  dye  introduced  to  the  system  through  the  center  port  was  also 
carefully  controlled  so  that  the  breakdown  form  was  clear  for  the  camera.  Excess  dye 
prohibited  accurate  measurements  of  the  stagnation  point,  as  did  not  enough  dye. 
Lighting  was  also  a  factor.  Medium  diffuse  back  lighting  was  used  to  illuminate  the 
vortex  core  and  breakdown  region. 

At  30  frames  per  second  it  was  nearly  impossible  to  separate  any  part  of  the 
structure  of  the  breakdown  from  the  turbulent  cone.  It  was  also  apparent  that  the  entire 
breakdown  structure  moved  as  the  stagnation  point  fluctuated,  such  that  there  was  no 
obvious  lag  time  or  remnants  from  the  previous  breakdown  locations.  This  observation  is 
merely  qualitative,  but  it  illustrates  that  when  the  CVB  occurs  over  a  delta  wing  aircraft, 
the  lift  and  drag  characteristics  of  the  aircraft  would  change  instantaneously  as  the 
breakdown  point  changed  location  on  the  wing.  To  further  refine  the  measurement  of  the 
CVB  stagnation  point  some  electronic  means  could  have  been  used,  but  such  an  effort 
was  not  assumed  because  visual  observation  was  possible.  A  potential  topic  for  a  future 
thesis  could  be  the  development  of  an  electronic  means  for  measuring  the  location  of  the 
stagnation  point. 
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C.  DATA  COLLECTION  AND  ANALYSIS 

The  video  tapes  were  played  back  using  a  Panasonic  AG- 1980  SVHS  Desktop 
Editor  and  a  Panasonic  TR-124MA  video  monitor  and  stopped  every  5  frames  to  record 
the  location  of  the  stagnation  point  of  the  CVB  using  a  ruler  mounted  to  the  monitor. 
2000  data  points  were  collected  in  this  manner  at  each  Reynolds  number.  Invariably  the 
axial  location  of  the  stagnation  point  of  the  breakdown  showed  at  least  some  movement 
in  between  each  of  these  measurements. 

It  is  necessary  to  note  that  there  were  several  breaks  in  the  video  which  prevented 
the  collection  of  2000  consecutive  data  points.  However,  since  it  seems  that  the 
fluctuations  are  at  all  times  random  this  should  not  affect  any  statistical  results.  It  is  also 
necessary  to  discuss  the  uncertainty  in  the  data  collection  process.  Because  the 
stagnation  point  of  the  vortex  breakdown  is  the  location  where  the  conical  form  erupts 
from  the  vortex  core,  the  exact  location  of  the  point  of  eruption  is  indetenninable  and  can 
only  be  estimated.  This  estimation  became  more  difficult  when  the  vortex  core  was 
enlarged  due  to  the  amount  of  dye  being  injected.  But  because  the  breakdown  was 
conical,  it  was  possible  to  trace  the  periphery  of  the  breakdown  toward  the  stagnation 
point  to  obtain  a  better  estimate  of  its  location.  The  uncertainty  in  these  measurements 
was  determined  to  be  ~2.0  mm  on  the  measurement  scale. 

The  mean  of  the  measured  fluctuations  was  calculated  and  then  subtracted  from 
each  data  point  to  give  the  fluctuations  from  the  mean.  Since  the  ruler  scale  was 
arbitrarily  placed,  it  was  necessary  to  convert  the  measured  fluctuations  from  the  ruler 
scale  to  actual  distances  for  comparison  with  future  studies.  Consequently  a  scaling 
factor  was  calculated  by  taking  the  ratio  of  a  distance  measured  on  the  ruler  to  the  same 
distance  measured  on  the  scale  in  the  video.  Using  the  scaling  factor,  each  measured 
fluctuation  was  subsequently  converted  to  actual  displacements  from  the  mean.  These 
fluctuations  of  the  CVB,  the  darting  of  the  stagnation  point,  were  finally  normalized  by 
the  tube  entrance  diameter.  Statistical  analysis  of  the  nonnalized  data  was  perfonned 
using  Microsoft  Excel  and  MATLAB  to  determine  the  physical  behavior  of  the 
phenomena.  The  results  are  presented  in  the  next  chapter. 
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III.  STATISTICAL  RESULTS 


A.  GENERAL  REMARKS 

This  chapter  deals  with  the  results  of  the  darting  of  the  stagnation  point  of  CVB  in 
tubes  and  the  distribution  of  these  fluctuations  for  comparison  with  investigations  carried 
out  with  delta  wings  and  subsequently  to  discuss  any  similarities  and  differences  of  the 
delta  wing  case  and  the  vortex  tube  case.  One  can  anticipate  that  the  position  of  the 
breakdown  as  well  as  the  distribution  of  the  fluctuations  in  the  two  said  cases  will  not 
necessarily  be  the  same.  The  fundamental  differences  are  that  in  the  case  of  swirl 
confined  to  a  tube,  total  circulation  throughout  the  length  of  the  tube  remains  constant, 
where  as  in  the  case  of  delta  wings  the  circulation  of  the  vortex  increases  along  its  length 
from  the  apex  to  the  tail  section  or  until  the  point  where  the  VB  occurs.  Thus  the 
differences  of  the  histories  of  the  circulation  are  expected  to  influence  not  only  the  shape 
of  the  distribution  of  the  stagnation  point  location  but  also  the  distance  at  which  the 
breakdown  occurs.  In  addition  to  the  foregoing,  one  must  point  out  additional  factors 
which  affect  all  of  the  characteristics  (breakdown  location,  frequency  of  fluctuation,  the 
distribution  of  the  fluctuations,  etc.)  in  the  two  cases. 

For  a  delta  wing  the  rolled  up  vortex  is  bounded  by  the  upper  surface  of  the  wing 
and  is  otherwise  free  in  all  directions;  the  external  pressure  gradient  remains  the 
atmospheric  pressure.  However  on  the  wing  surface  the  swirling  conical  flow  gives  rise 
to  oppositely-signed  vorticity.  Therefore  the  vorticity  generated  on  the  lower  surface  of 
the  wing  and  curling  up  to  a  vortex  on  the  upper  surface  increases,  as  stated  earlier,  with 
distance,  but  is  also  decreased  (to  a  smaller  extent)  due  to  the  vorticity  generated  at  the 
upper  surface.  After  nearly  50  years  of  research  on  VB  it  has  not  been  possible  to 
analytically  predict  the  position  of  the  breakdown,  let  alone  its  consequences.  One  can 
certainly  state  that  the  flow  over  the  delta  wing  can  not  be  the  same  as  in  a  tube.  From  a 
practical  point  of  view  a  delta  wing  is  seldom  free  from  protuberances,  particularly  at  its 
upper  surface,  which  interfere  with  the  conical  flow. 

In  the  case  of  a  tube,  the  circulation  is  constant  even  though  the  vorticity 
distribution  varies  first  gradually  as  one  approaches  the  stagnation  point  and  rapidly  at 
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and  beyond  the  stagnation  point  as  evidenced  in  the  photographs  shown  in  Figure  5.  The 
tube  walls  are  sufficiently  far  away  from  the  vortex  core  (prior  to  breakdown)  to  account 
for  the  changes  in  vorticity  and  the  effect  of  the  pressure  gradient  is  judiciously 
compensated  by  the  slight  expansion  of  the  tube  wall.  Thus  this  makes  the  investigation 
of  the  fluctuations  of  VB  in  a  tube  relatively  ideal  and  more  physics  based  without  trying 
to  deal  with  other  influencing  factors  such  as  the  changes  in  the  vorticity  distribution,  etc. 
as  noted  above.  Consequently  the  creation  of  a  mobile  stagnation  point  and  its 
fluctuation  created  by  the  occurrence  of  the  breakdown,  as  better  understood  by  the 
distribution,  cannot  be  predicted  because  the  transformation  of  a  vortex  with  a  very  thin 
core  into  a  swirling  turbulent  vortex  with  a  large  core  and  highly  distributed  vorticity 
with  numerous  orders  of  coherent  structures  leads  to  not  only  one  but  to  a  series  of 
discrete  frequencies. 

Comparison  of  the  results  of  the  delta  wing  studies  with  those  obtained  in  a  tube 
are  far  from  obvious  and  had  to  be  understood  towards  the  understanding  of  the 
occurrence  of  the  breakdown.  One  may  conjecture,  for  example,  that  the  breakdown 
location  will  exhibit  both  the  high  frequency  fluctuations  and  then  several  lower 
frequencies  due  to  the  response  of  the  large  turbulent  fluctuations  of  flow  down  stream  of 
the  breakdown.  In  fact,  cold-flow  experiments  conducted  in  combustion  chambers  have 
shown  that  the  shorter  the  length  of  the  combustion  chamber  for  a  given  diameter  the 
larger  the  range  of  frequencies  (Swithenbank  and  Chigier  1969).  In  the  present 
experiments  the  full  awareness  of  the  numerous  phenomena  noted  above  dictated  that  the 
length  of  the  tube  downstream  of  the  breakdown  point  be  made  sufficiently  long  (L/D  ~ 
27)  to  eliminate  Helmholtz  resonance  and  additional  frequencies.  Furthermore  the  exit  of 
the  tube  was  varied  numerous  times  to  reduce  such  fluctuations.  An  additional 
consideration  for  the  downstream  section  of  the  tube  was  that  the  conical  breakdown  fill 
the  entire  tube  long  before  the  end  of  the  tube  was  reached  so  that  essentially  the  flow 
became  self-similar.  It  is  on  the  basis  of  the  foregoing  that  we  will  now  discuss  the 
results  and  compare  them  with  those  from  delta  wings. 
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B.  DISCUSSION  OF  FINDINGS 

As  expected  the  external  influences  on  the  evolution  of  the  vortex  over  a  delta 
wing  exhibit  fluctuations  which  are  close  enough  to  a  Gaussian  distribution  but  also 
exhibiting  wilder  fluctuations  along  the  length  of  the  vortex.  It  is  a  well  known  fact  that 
part  of  the  secondary  fluctuations  come  from  the  variation  of  the  breakdown  location 
with  the  tail  end  of  the  wing.  In  other  words  when  the  VB  momentarily  moves 
downstream  from  its  mean  position  it  comes  closer  to  the  tail  of  the  wing  and  thus  is 
further  influenced  by  its  proximity.  When  the  breakdown  moves  upstream  (for  random 
vorticity  fluctuations)  then  the  increase  of  the  distance  between  the  breakdown  and  the 
tail  reduces  the  larger  fluctuations.  Evidently  the  extraneous  influences  are  almost 
entirely  eliminated  in  the  case  of  our  experiments.  For  example  Menke’s  (1996) 
probability  density  function  of  fluctuations  of  breakdown  locations  over  delta  wings  for 
angles  of  attack  of  26  and  37  degrees  are  shown  in  Figure  10(a)  and  (b)  respectively. 


(a) 
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(b) 

Figure  10.  Probability  Density  Function  of  fluctuations  of  breakdown  for:  (a)  26 
degree  and  (b)  37  degree  angle  of  attack.  (After  Menke  1996) 

The  figure  shows  that  the  data  is  skewed  toward  the  right  of  the  Gaussian  distribution 
(shown  with  dotted  lines)  at  the  relatively  smaller  angles  of  attack.  However  at  very 
large  angles  of  attack  the  tail  proximity  effects  are  so  strong  that  the  distribution  comes 
not  so  surprisingly  closer  to  a  Gaussian  distribution  with  at  least  three  wild  fluctuations. 

In  some  cases  the  fluctuation  distribution  is  not  presented  as  a  probability  density 
function,  but  rather  in  the  fonn  histograms  of  the  dominant  frequencies,  the  power 
spectra  density.  Figure  11  shows  Mitchell’s  (2000)  presentation  of  his  power  spectra 
density  results  as  compared  to  those  of  Menke  (1996). 
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I _ L 

10  s 

Strouhal  Number 

Figure  11.  Fluctuations  of  breakdown  location:  Mitchell’s  results  are  labeled 

“current  results”  and  are  superimposed  on  Menke’s  (1996)  results.  (After 

Mitchell  2000) 

The  overlap  of  the  results  in  these  delta  wing  experiments  are  for  Strouhal  numbers  of 
about  0.04  to  0.1.  Gursul  (1995)  also  has  comparable  results  with  peaks  in  the  power 
spectra  density  at  Strouhal  numbers  of  0.1,  0.06,  and  0.08  for  angles  of  attack  of  37,  30 
and  27  degrees  respectively.  When  compared  to  the  results  from  the  present  study  in 
Figure  7  the  correlation  between  the  dominant  frequency  of  fluctuations  in  a  tube  and 
over  delta  wings  is  evident  for  the  stated  Strouhal  number  range. 

We  will  now  discuss  the  results  obtained  in  the  present  investigation  in  light  of 
the  foregoing.  A  critical  investigation  of  the  valuable  statistical  properties  of  CVB  must 
start  with  an  understanding  of  the  time  history  of  the  breakdown  location.  Observation  of 
video  of  the  CVB  reveals  that  the  stagnation  point  of  the  breakdown  is  always  changing 
location.  As  shown  in  Figure  12,  the  fluctuation  has  no  apparent  pattern  and  moves 
arbitrarily  about  the  mean  value,  in  this  case  the  x-axis. 
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Figure  12.  Time  history  of  breakdown  locations  for  ReD=230,000.  500  data  points 

(Top)  and  1,000  data  points  (Bottom) 


Figure  13  shows  the  results  for  the  normalized  fluctuation  from  the  mean  as 
compared  with  the  Gaussian  distribution. 


(a) 
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1.0 


.  ReD=300,000  Data  Distribution 


(b) 

Figure  13.  Comparison  of  normalized  data  to  a  Gaussian  distribution  (a) 

ReD=23 0,000  (b)  ReD=300,000. 

Evidently  these  distributions  are  almost  perfectly  Gaussian  except  at  the  very  tail 
ends  showing  the  larger  fluctuations  are  more  apt  to  contain  other  frequencies  than  the 
smaller  fluctuations.  Figure  13  shows  the  non-dimensional  fluctuation  from  the  mean,  p, 
where  p=5.7 1  e-4  for  ReD=230,000  and  p=8.63c- 1 6  for  ReD=300,000.  The  standard 
deviation,  a,  was  a=0.363  for  Reo=230,000  and  a=0.451  for  ReD=300,000.  Figure  14 
shows  that  the  cumulative  probability  distribution  for  both  Reynolds  numbers  is  very 
symmetric. 
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Fluctuation  from  the  Mean:  X'/D 


(a) 


(b) 

Figure  14.  Cumulative  Probability  Distribution  for  normalized  fluctuations  from  the 

mean,  (a)  Reo=230,000  and  (b)  ReD=300,000. 


The  experiments  were,  as  noted,  carried  out  at  two  significantly  different  Reynolds 
numbers.  The  results  show  that  at  such  high  Reynolds  numbers  the  effect  of  the 
Reynolds  number  is  insignificant  when  protuberances  are  absent. 
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For  a  more  complete  understanding  of  the  data  presented  here,  the  fluctuations 
will  now  be  presented  in  the  form  (Xm+X')/D  or  X/D,  where  Xm  is  the  mean  stagnation 
point  location,  X'  is  the  fluctuation  from  the  mean,  and  D  is  the  entrance  diameter  of  the 
tube  as  shown  in  Figure  15. 


Figure  15.  Tube  test  section  indicating  the  location  of  the  mean  stagnation  point,  Xm, 

and  the  entrance  diameter,  D. 


The  distribution  of  the  data  for  both  Reynolds  numbers  is  shown  in  Figure  16. 


Fluctuation  from  the  Mean:  X/D 

(a) 
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Fluctionation  from  the  Mean:  XI D 

(b) 

Figure  16.  Comparison  of  normalized  data,  X/D,  to  a  Gaussian  distribution  (a) 

ReD=230,000  (b)  ReD=300,000. 

Note  that  the  mean  and  standard  deviation  change  to  p=7. 1 34  and  <7=0.556  for 
ReD=230,000  and  ju=7.176  and  <7=0.623  for  ReD=300,000.  The  shapes  of  the 
distributions  are  nearly  Gaussian  and  both  are  slightly  skewed  to  the  left.  These  results 
do  however  resemble  the  previously  presented  results  from  Figure  13  and  support  the 
conclusion  that  the  darting  of  the  stagnation  point  is  normally  distributed. 

C.  CONCLUSIONS  AND  RECOMMENDATIONS 

It  has  been  shown  that  conical  vortex  breakdown  exhibits  streamwise  fluctuations 
in  the  location  of  the  stagnation  point  that  are  arbitrary  and  naturally  occurring.  The 
frequencies  at  which  these  fluctuations  occur  in  tubes  are  identical  to  those  found  in  the 
delta  wing  cases,  validating  this  study  as  a  positive  contribution  to  the  bank  of  knowledge 
to  include  other  such  phenomena  and  physical  attributes  of  vortex  breakdown.  After 
nearly  50  years  of  research  on  vortex  breakdown,  to  include  theoretical,  numerical,  and 
experimental  investigations,  it  is  the  case  that  this  phenomenon  has  been  limited  to  and 
remains  largely  in  the  realm  of  qualitative  observation.  The  complex  geometries  and 
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protuberances  in  delta  wing  aircraft  make  it  difficult  to  accurately  predict  the 
instantaneous  location  of  the  stagnation  point  and  since  the  history  of  the  breakdown 
plays  a  decisive  role  in  the  movement  of  the  breakdown  as  well  as  the  interaction  of  the 
two  leading  edge  vortices,  which  are  formed  on  the  wings  of  the  delta  wing,  prediction  of 
the  breakdown  location  is  at  this  time  too  complex  and  involves  too  many  variables. 
Experiments  in  tubes  have  provided  robust  results  and  should  be  depended  upon  for 
comparison  with  future  numerical  simulations.  As  computational  methods  become  more 
advanced  and  models  are  developed  to  predict  high  Reynolds  number,  turbulent,  swirling 
flows  the  physical  limitations  in  measurement  techniques  of  vortex  breakdown  should  no 
longer  be  an  obstacle.  However,  until  this  definitive  point  is  reached,  the  efforts  of 
innumerable  researchers  using  various  methods  will  encompass  the  basis  of 
understanding  of  vortex  breakdown. 
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